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ABSTRACT 

We present optical integral field spectroscopy of the Ha-luminous (> 10^^ ergs~^ ) central 
cluster galaxies in the cores of the cooling flow clusters A1664, A1835, A2204 and Zw8193. 
From the [NII]+Ha complex in these moderate resolution (70-150 km ) spectra we derive 
two-dimensional views of the distribution and kinematics of the emission line gas, and further 
diagnostics from the [SII] and [01] lines. 

The Ha emission shows a variety of disturbed morphologies, ranging from smooth but 
distorted to clumpy and filamentary, with velocity gradients and splittings of several hun- 
dred kms~^ on spatial scales of 20kpc or more. Despite the small sample size, there are 
some generic features. The most disturbed Ha emission appears to be associated with sec- 
ondary galaxies within 10-20kpc (projected) of the central galaxy and close in velocity to the 
Ha. The global Ha kinematics match those of the CO(l-O) emission in single-dish data. The 
[NII]/Ha, [SII]/Ha and [OI]/Ha ratios vary little with position, local Ha surface brightness 
or between clusters. 

We propose that the Ha and CO emission arise in molecular clouds heated by a starburst, 
and that the latter has been triggered by interaction with a secondary galaxy. Such CO emis- 
sion is known to trace massive (> 10^^ M©) compact (< 20 kpc) reservoirs of cool molecular 
gas, and it is plausible that an inf ailing galaxy would disturb this gas, distorting the Ha mor- 
phology and initiating widespread star formation. We also examine the role of cloud-cloud 
collisions in the undisturbed molecular gas reservoir, and suggest that they might be an im- 
portant source of excitation for the emission line gas in the cores of lower Ha luminosity 
clusters with less intense star formation. 

Key words: galaxies: clusters: individual: A1664, A1835, A2204, ZW8193 - cooling flows - 
intergalactic medium - ionized gas 



1 INTRODUCTION 

In recent years our understanding of the X-ray cooling flow phe- 
nomenon in galaxy cluster cores has been revolutionised. Through- 
out most of the 1980s and 1990s, X-ray observations suggested that 
gas in the central 100 kpc is cooling out at rates of up to several hun- 
dred solar masses per year, but the lack of evidence for a reservoir 
of cooled gas led to heated debate (summarised by Fabian 1994) 
over this interpretation of the X-ray data. Results from XMM- 
Newton and Chandra have since led to a sharp downward revision 
in X-ray cooling rates (e.g. Schmidt, Allen & Fabian 2001) and 
also reveal a strong deficit of line emission from gas cooling be- 
low Tviriai/3 (Peterson et al. 2003). The implication is that X-ray 
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cooling is quenched, for which numerous mechanisms have been 
proposed, including: rapid mixing of hot and cold phases, inho- 
mogeneously distributed metals in the intracluster medium (Fabian 
et al. 2001,2002); active galactic nucleus (AGN) heating by jets 
(Briiggen & Kaiser 2003) and sound waves (Fabian et al. 2003); 
thermal conduction of heat from the hotter outer parts of the clus- 
ter into the cooler core (Voigt et al. 2002); a significant relativistic 
cosmic ray component frozen into the thermal gas (Cen 2005); the 
release of gravitational energy from blobs of gas which detach from 
the bulk flow and fall directly into the core (Fabian 2003). 

Concurrently, significant progress has been made in identi- 
fying cool gas and dust in cluster cores. Edge (2001) detected 
CO emission in the centres of 16 cooling flows, consistent with 
10^ - 10^^-^ M© of Hs at 20-40 K for a standard C0:H2 con- 
version (see also Salome & Combes 2003). These are roughly the 
masses expected, given the revised cooling rates and likely ages. In- 
terferometry shows further that the CO emission is localised within 
the central few arcsec of the cluster (Edge & Frayer 2003; Sa- 
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lome & Combes 2004). The frequent occurrence of smaller masses 
(~ 10^ - 10^ Mq) of hot H2 has also been estabhshed (e.g. Edge 
et al. 2002; Jaffe, Bremer & van der Werf 2001), and excitation 
analysis suggests that this hot H2 is a high pressure, transiently- 
heated component (Wilman et al. 2002). Both CO and H2 emis- 
sions correlate well with the strength of the Ha emission from ion- 
ized gas at 10^ K, whose prevalence in these environments, often 
in the form of spectacular filaments, has long been known (e.g. Hu 
et al. 1983; Crawford et al. 1999). Despite the clear association be- 
tween optical line emission and short central X-ray cooling times 
(Peres et al. 1998; Bauer et al. 2005), their physical relationship 
is ill-understood. Photoionisation by the observed excess popula- 
tion of hot massive stars can energetically account for the Ha lu- 
minosities in the most luminous systems (Allen 1995; Crawford 
et al. 1999). Spitzer MIPS photometry of 11 CCGs by Egami et 
al. (2006) also shows that the most Ha-luminous in their sample 
(A 1835, A2390 and Zw3146) have prominent far-infrared thermal 
dust emission plausibly powered by star formation, two of them 
with LiR > 10^^ L0. At lower Ha luminosities the picture is less 
clear: the tapping of energy from the intracluster medium (ICM) 
through turbulence (Crawford & Fabian 1992) and heat (Sparks et 
al. 2004) are just two mechanisms which have been invoked to ex- 
plain the optical nebulosity in such systems. 

In this paper we present integral field spectroscopy of the 
ionized gas in the cores of four such clusters, A1664, A1835, 
A2204 and Zw8193. The principal aim is to obtain a full two 
dimensional view of the distribution and kinematics of the gas 
through high resolution spectroscopy of the Ha-i-[NII] emission 
line, with additional ionization information being gleaned from 
the [SII]AA6717,6731 and [OI]AA6300,6363 lines where possible. 
These four central cluster galaxies (CCGs) all have Ha luminosi- 
ties exceeding 10^^ ergs~^ , making them 4 of the top 6 most- 
Ha luminous systems in the extensive CCG spectroscopic survey 
by Crawford et al. (1999). In this regime of Ha luminosity, pho- 
toionisation by a young stellar population can account energetically 
for the luminosity of the Ha nebulosity (Allen 1995; Crawford et 
al. 1999). In addition to an analysis of the CCGs, we also present 
spectroscopy of other sources within the IFU field of view, includ- 
ing other cluster galaxies and (in the case of A2204) a serendipitous 
gravitationally-lensed background galaxy. We first present results 
for the clusters individually and then summarise and interpret their 
generic features. Throughout the paper we assume a cosmology 
with Ho = 70kms"^ Mpc"\ Qm = 0.3 and Qa = 0.7 and all 
physical quantities quoted from other papers have been converted 
accordingly. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 VLT-VIMOS Integral Field Spectroscopy 

The observations of A 1664, A 1835 and A2204 were taken in ser- 
vice mode on 2003 April 1 1 with the integral field unit (IFU) of the 
Visible Multiobject Spectrograph (VIMOS) on UT3 of the 8.2m 
Very Large Telescope (VUT) at ESO Paranal in Chile (for further 
information on VIMOS see LeFevre et al. 2003). They were among 
the first taken with the VIMOS IFU in open time. The IFU was op- 
erated in HR-Red mode, offering a 27'' x 27'' field of view covered 
by 1600 optical fibres of 0.67" diameter. The fibres are coupled to 
a microlenses to ensure near-continuous sky coverage. The field of 
view in this IFU mode is split into 4 quadrants of 400 fibres, three of 
which disperse the light with the HR-Red grism over 6250-8700A, 



while the fourth quadrant employs the HR-Orange grism spanning 
5500-7450A. The dispersion and spectral resolution are approxi- 
mately 0.6A per pixel and 1.8 A FWHM, respectively. For each of 
the three targets a pair of 1670s exposures was taken, with a point- 
ing dither of (Aa, A^) = (—6", —6") between them. The seeing 
was in the range 0.5-1" throughout. Table 1 provides an observa- 
tion log. 

The data were reduced on site at the National Institute for As- 
trophysics (INAF) in Milan using the VIMOS Interactive Pipeline 
and Graphical Interface (VIPGI) (Scodeggio et al. 2005) developed 
by the VIRMOS consortium. The raw two-dimensional frames for 
each quadrant were bias-subtracted, and the spectra for the indi- 
vidual fibres were traced, extracted and flat-fielded. The pipeline 
wavelength calibration was affected by instrumental flexure in the 
time interval between science and calibration observations. The 
wavelength calibration was thus performed manually in IRAF us- 
ing night-sky emission features, rebinning to a common dispersion 
of 0.6 A per pixel. Using the mapping from detector to aperture fibre 
position, the four quadrants were united into a single 40 x 40 fibre 
data cube for each pointing, covering the wavelength range 5500- 
8740A (i.e. large enough to accommodate the wavelength range of 
the HR-red and HR-orange grisms). The raw data format is such 
that the centres of adjacent fibres are separated by 5 pixels on the 
detector whilst the light from each fibre covers 3 pixels FWHM. 
There is thus a small amount of fibre-to-fibre light leakage imply- 
ing that the spectra of individual fibres will not be completely inde- 
pendent. However, since the seeing is comparable to the fibre size 
and we are chiefly concerned with diffuse extended emission, the 
effect is not considered important and no corrections for it were 
made. 

Subtraction of the night- sky emission background was per- 
formed by fitting gaussian profiles to a triplet of sky features at 
7316, 7329 and 7341 A in each fibre. This intensity information was 
used to correct for fibre-to-fibre transmission variations across the 
field of view, and then to create a night-sky emission spectrum from 
a source-free region which was then subtracted across the whole 
cube. Due to the use of different grisms, sky subtraction was per- 
formed separately for quadrants 1-3 and 4, respectively. For A1664 
and A2204 the two dithered sky- subtracted datacubes for each tar- 
get were then combined into a single datacube 49 x 49 fibres in size 
(33" X 33"). For A1835, only the second of the two pointings con- 
tains significant line emission, since the target was inadvertently 
placed too near the edge of the field of view. Flux calibration of the 
spectra was not performed, but comparison with published emis- 
sion line slit fluxes (Crawford et al. 1999) suggets that Ha emission 
is detected down to minimum surface brightness levels in the range 
0.5-1.0 X 10~^^ erg cm~^ s~^arcsec~^ . Following the reduction 
of the raw data, all subsequent cube manipulation and emission line 
analysis were performed within IDL. 

2.2 WHT-OASIS Integral Field Spectroscopy 

Observations of the redshifted Ha emisson in Zw8193 were made 
with the OASIS IFU on the William Herschel Telescope (WHT) on 
La Palma on 2005 September 2 in 0.8" seeing and photometric con- 
ditions. We used the MR 807 filter which covers a wavelength range 
of 7690 - 8460A at a spectral resolution of A/AA=2020 (2.2A per 
pixel). The 7.4" x 10.3" field of view was sampled contigiously 
with 0.26" fibres. The total integration time was 7200s, split into 
4x30mins and each dithered by one IFU lenslet to account for bad 
pixels. For further information about the OASIS IFU at the WHT 
see Benn et al. (2003). 
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Table 1. Summary of the observations 



Target 


z 


IFU 


Exposure 


Seeing 


Scale 








time (s) 


(arcsec) 


kpc/arcsec 


A1664 


0.1276 


VIMOS 


3340 


0.5-1.0 


2.3 


A1835 


0.2523 


VIMOS 


3340 


0.5-1.0 


3.9 


A2204 


0.1514 


VIMOS 


3340 


0.5-1.0 


2.6 


Zw8193 


0.1825 


OASIS 


7200 


0.8 


3.1 



The OASIS data were reduced using the XOASIS data- 
reduction package (Rous set 1992) developed at CRAL (Lyon). The 
reduction steps include bias and dark subtraction, extraction of the 
spectra using an instrument model, wavelength calibration, cosmic- 
ray rejection, low frequency flatfielding and sky- subtraction. The 
individual data-cubes were then mosaiced (taking into account the 
small spatial offsets between exposures) and resampled onto a com- 
mon grid of 0.26'^ x 0.26^' to create the final datacube. The limit- 
ing emission line surface brightness is comparable to that reached 
in the VIMOS observations. Line fitting to the OASIS data was 
performed by binning up 3 x 3 pixels in order to better match the 
seeing. 



3 A1664 

3.1 Ha morphology and kinematics 

With a extinction-corrected Ha slit luminosity of 1.6 x 
lO^^ergs"^ , the central cluster galaxy of A 1664 is among the 
brightest Ha emitters in the optical spectroscopic survey of the 
ROSAT Brightest Cluster Sample (BCS; Crawford et al. 1999); fits 
to the excess blue continuum light imply a visible star formation 
rate of 23 M© yr"^ (see also Allen 1995). 

The complexity of the line emission in A1664 is evident 
in the reduced two-dimensional frame of one of the quadrants, 
shown in Figure 1. The emission spans a large range in veloc- 
ity and at many spatial locations two velocity components are re- 
quired to satisfactorily model the Ha-h[NII]AA6548,6584 profiles. 
Within each spatial pixel, the lines were fitted with zero, one or 
two Ha-i-[NII]AA6548,6584 gaussian velocity components atop a 
flat continuum level; within each velocity component the Ha/[NII] 
was allowed to vary but the [NII]AA6548,6584 double ratio was 
fixed at the theoretical 1:3 ratio. A simple chi-squared fitting cri- 
terion was used to decide whether one, two or zero velocity com- 
ponents were fitted at each position. The complex was fitted within 
the wavelength range 7330-7450A (i.e. 200 dispersion bins) and 
the noise-level was estimated from the spectrum in the wavelength 
range 6490-7500A. Chi-squared values were evaluated for fits with 
0, 1 or 2 velocity components, Xo^ Xi ^^d xi. respectively. For 
Xo~xi > 150, 1 velocity component was fitted; for Xi~X2 > 10^ 
two components were deemed necessary. All fits were examined by 
eye to confirm whether they were reasonable and to remove spuri- 
ous fits to obvious noise spikes. 

Fig. 121 shows the Ha and continuum maps derived from these 
emission line fits. For clarity, three separate Ha maps are shown: 
the total flux (i.e. all the Ha emission), and the separate red and 
blue component fluxes for those pixels fitted with two Ha velocity 
components. The velocity and FWHM fields derived from these fits 
are shown in Fig.|3] and for clarity are again shown separately for 
the red and blue components where two components are required, 
and for those regions where only a single velocity component is 



needed. Line profile fits are shown in Fig.Qfor selected regions of 
Ha emission in Fig.|2| From Fig.|2| we see that the Ha emission has 
a linear filamentary structure spanning approximately 13 arcsec (31 
kpc) and passing through the nucleus of the central cluster galaxy 
along a north-east - south-west axis. Two velocity components are 
detected along its whole length: the red component peaks on the 
nucleus of the CCG, but the blue component exhibits two maxima 
(labelled Bl and B2) located approximately 1.5 arcsec from the 
nucleus, along an axis perpendicular to the general direction of the 
filament. This filament and nuclear bipolar structure are discussed 
further below. 



3.2 The 30-kpc Ha filament and nuclear bipolar structure 

The presence of two velocity components along the length of the 
filament (most clearly at its NE end where complications from the 
nuclear emission of the CCG are minimal) suggests that it might be 
due to a disturbed secondary emission component superimposed 
on a quiescent background, or a stream of gas wrapped around the 
CCG, or possibly related to some form of outflow. It is plausibly re- 
lated to the galaxy (obj #1) situated just beyond the SW extremity 
of the filament in Fig.|2| which is confirmed to be at the cluster red- 
shift (see Fig.|3} with a line of sight velocity of —560 km s~ ^ in the 
velocity system adopted in Fig.|3] This gas may have been stripped 
from or disturbed by obj #1. Indeed, the velocity of the blue emis- 
sion component at position B2 (—450 km ) is similar to that of 
obj#l (-560 km s-^- 

The bipolar structure may instead be an outflow from the 
CCG driven by starburst activity. The line of sight velocity dif- 
ference between the blue emission components at Bl and B2 is 
~ 350 kms~^ . This suggests a minimum deprojected outflow ve- 
locity of 175 kms~^ , but the symmetrical location of Bl and B2 
about the nucleus suggests that the outflow may be principally di- 
rected into the plane of the sky with a much higher velocity. In the 
latter case, the observed radial extent of the outflow (1.5^^=3.5 kpc) 
can be used to constrain its age. Assuming a deprojected outflow 
speed of 600kms~^ as for the Ha filaments in M82 (Shopbell 
& Bland-Hawthorn 1998), the implied age is ~ 6 Myr. By con- 
trast, the orbital timescale of galaxy #1 responsible for the creation 
of the filament is ~ 100 Myr (assuming a circular orbit of ra- 
dius 15 kpc traversed at 500kms~^ ). Considering the mis-match 
in these timescales, it is too coincidental that an outflow from a 
nuclear starburst triggered by the passage of galaxy #1 should be 
observed to have a similar extent to the overall stream of gas. We 
thus favour an interpretation in which the complex Ha structure is 
due to a disturbance induced by, or the stripping of gas from, galaxy 
#1. 



3.3 Archival VLT-FORSl V, R and I imaging 

To shed further light on the origin of the complex Ha distribution, 
we retrieved archival VLT-FORSl imaging of A 1664. The dataset 
comprises images three 330 second exposures taken through the 
V, R, and I filters (Bessel photometric system) on 2001 May 29 in 
photometric conditions and 0.8^' seeing. Following bias- subtraction 
and flat-fielding, we merged them into a true-colour image which 
is displayed in Fig.|9l Ha contours from the VIMOS IFU observa- 
tions were aligned and overlaid on this image. 

The main features apparent from Fig.|9lare that: (i) the north- 
eastern quadrant of the CCG is redder than other parts of the CCG, 
suggesting dust obscuration; (ii) the south-western quadrant of the 
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CCG is bluer than other cluster galaxies, suggesting more recent 
star formation; (iii) the red Ha velocity component peaks on a nu- 
clear dust lane which bisects the two central peaks of the blue Ha 
velocity component; (iv) the southern extension of the Ha emis- 
sion (labelled B4 in Fig.|2j coincides with a blue knot of continuum 
emission. 

V-R and R-I colours (expressed relative to the colours of the 
nearest A 1664 cluster elliptical) were extracted for regions B 1,2,4 
and object #1 in Fig.|2| These show that (i) galaxy #1 has the same 
colour as other cluster ellipticals at this redshift; (ii) the CCG is 
~ 0.1 mag bluer m V — R and ~ 0.1 mag redder m R — I 
and that the intrinsic dust extinction increases across it from south- 
west to north-east. From long- slit spectra, the intrinsic extinction 
inferred from the B aimer decrement is in the range 0.46-0.63 mag 
(Crawford et al. 1999; Allen 1995). Corrected for this extinction, 
the same authors show that approximately 90 per cent of the (rest- 
frame) 4500A continuum is due to 05 and B5 stars giving rise to 
the 23 M0 yr~^ starburst. The gradient of increasing dust extinc- 
tion from sout-west to north-east could be due to the gradient in 
time elapsed since the passage of galaxy #1, which may have trig- 
gered star formation and dust production. U or B-band imaging 
would more clearly show the distribution of recent star formation. 

3.4 Ionization properties 

Although the VIMOS data contain the emission lines of 
[SII]AA6717,6731 and [OI]AA6300,6363 in addition to [NII]-hHa, 
we do not present an exhaustive study of the ionization properties 
of the gas. Instead, we impose the kinematics found for [NII]-i-Ha 
and accordingly fit the [SIX] and [01] doublets with zero, one or 
two velocity components with only the line intensities allowed to 
vary, i.e. we do not search for emission in these lines where there is 
no [NII]-i-Ha detection. In so doing, we discover that the ionisation 
properties of the Ha gas are quite uniform, with no notable spa- 
tial gradients. In Fig.|5]we show scatter plots of line ratios versus 
Ha flux for the combined sample of lenslets. There is essentially no 
variation in [Nil] A6584/Ha and [01] A6300/Ha over more than one 
decade in Ha surface brightness. These line ratios are consistent 
with the integrated slit values found by Crawford et al. (1999). The 
uniformity of the ionization state across the source supports our 
interpretation that all the gas (the filament and the B1/B2 nuclear 
structure) has the same origin. For comparison, the emission line 
filaments in the kiloparsec- scale superwind outflow in M82 have 
[NII]A6584/Ha = 0.3 - 0.6 (Shopbell & Bland-Hawthorn 1998). 

3.5 Comparison with CO and H2 observations 

In the CO line survey of cooling flow CCGs by Edge (2001), 
A 1664 exhibits unusually broad CO(l-O) emission spanning 
620kms~^ full- width at zero intensity, with a total inferred mass 
of cool molecular hydrogen is 4.4 it 0.7 x 10^° M0. Consider- 
ing the 23.5 arcsec beam size for these IRAM 30-m observations, 
its kinematics are a close match to those of Ha; in particular, the 
skewness of the CO line towards negative velocities suggests that 
the blue Ha velocity components may have associated CO. A direct 
comparison is shown in Fig.|8] the CO(l-O) profile is taken from 
Edge (2001), and the Ha profile is formed by integrating over the 
whole source and adding synthetic gaussians for the Ha emission 
only, using the parameters derived from the multi-component fits 
to Ha-i-[NII]. This procedure is used in order to compare the kine- 
matics without complications from the blended [Nil] lines. Fig.|8] 
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Figure 4. The spectrum of continuum object # 1 in Fig.|2| The absorption 
feature due to the Na D A5892.5 blend is labelled, from which a redshift of 
z = 0.1255 is determined. The feature at 686O-688OA is telluric absorp- 
tion. 

demonstrates that the widths and centroids of the Ha and C0(1- 
0) emission are indeed comparable, and that much of the breadth 
of these lines can be attributed to the regions of kinematically dis- 
turbed Ha emission with two velocity components. 

For comparison, in the K-band UKIRT CGS4 spectrum 
of Edge et al. (2002), the width of the Paa line is 
1100 kms"^ FWHM, for a N-S slit of 5.5 x 1.2 arcsec. This is 
consistent with the observed Ha kinematics when the VIMOS data 
are extracted over the same area and fitted with a single velocity 
component. In contrast, the H2 v=l-0 S(3) line is not resolved, sug- 
gesting an intrinsic width much below the instrumental broadening 
of 570 km s"^ FWHM. This suggests that the hot H2 is distinct 
from the CO and Ha-emitting gas in the case of A1664. 



3.6 Cluster X-ray emission 

A1664 has been observed with the ACIS-S instrument on the Chan- 
dra X-ray observatory for 10 ks. In Fig.|6|we show the smoothed 
0.5-5 keV image (courtesy of A.C. Fabian). This shows that the X- 
ray emission is also elongated along a NE-SW axis with a possible 
filamentary structure, resembling the Ha emission but on a much 
larger spatial scale. A full analysis of the X-ray dataset will be pre- 
sented elsewhere. The observed similarity in Ha and X-ray mor- 
phology is analogous to that seen in the core of A 1795 (Crawford, 
Fabian & Sanders 2005), although in that cluster the emission line 
gas in the filament is very quiescent, in stark contrast to the case of 
A 1664. As Fig.|6|demonstrates, the CCG's small neighbour galaxy 
(labelled object #1 in Fig.|2| appears to coincide with a local mini- 
mum in the X-ray emission. 



4 A1835 

At a redshift of z = 0.2523, A1835 is the most distant cluster pre- 
sented in this study. It is also our most luminous Ha emitter, with 
an extinction-corrected luminosity of 3 x 10^^ ergs~^ and a vis- 
ible star formation rate of ^ 100 M© yr~^ (Crawford et al. 1999; 
Allen 1995). The Spitzer MIPS photometry of Egami et al. (2006) 
reveal far infrared thermal dust emission with Z/ir — 7 x 10 L©, 
plausibly powered by star formation. In the X-rays, the cluster is 
the most luminous in the ROSAT BCS and subsequent Chandra 
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Figure 1. Part of the reduced 2-D frame for a single quadrant of the A 1664 observation. The complexity and large velocity range of the 
emission in Ha+[NII]AA6548,6584 and [SII]AA6717,6731 are clearly evident. Note that neighbouring fibres on the sky do not necessarily 
lie in adjacent rows of this frame. 
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Figure 2. Ho; and continuum (6700-7000A observed frame) intensity plots from multi-component fits to the Ho; -i- [Nil] complex in A 1664. These images have 
been smoothed with a gaussian of FWHM=1.5 pixels (=1 arcsec). The emission line complexes were fitted with 1 or 2 velocity components, shown as grey 
and black, respectively. The Ha (BLUE/RED) panels show the intensities of the two components where fitted. Ha (TOTAL) shows the total Ha intensity and 
includes those pixels in which only one velocity component was fitted. The object labelled #1 in the continuum image is a galaxy at redshift z = 0.1255 (i.e. 
in the cluster core), as shown by its spectrum in Fig.|4| Line profile fits to regions B1-B4 are shown in Fig.0 
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Figure 3. Ha kinematics as derived from multi-component fits to the Hq;-i-[NII] complex in A1664. All units are kms~^ and the velocity fields are with 
respect to a common zeropoint of the nominal cluster redshift {z = 0.1276). The velocity range spanned by the colour bar varies from panel to panel in order 
to clearly show the kinematic features of interest. For clarity, the velocity and FWHM fields are shown separately for the red and blue velocity components 
where two kinematic components were required to fit the line, and for those fibres where just a single velocity component was required to fit it. For comparison, 
the velocity of continuum object #1 in Fig.|2|in this reference system is -560 km (determined from spectrum in Fig.|4|. For ease of comparison, the total 
Ho; intensity map of Fig.|2|is also shown. 



data are consistent with a young (6 x 10 yr) cooling flow operat- 
ing within r = 30 kpc (Schmidt, Allen & Fabian 2001). 

In comparison with A 1664 the kinematics are simple and rela- 
tively quiescent, requiring a single velocity component at each po- 
sition. The Ha map and velocity field are shown in Fig. 1101 The Ha 
emission has a maximum extent of 8^^ (31 kpc) and is elongated 
in a NW-SE direction; the reconstructed continuum image shows 
a small companion galaxy 5^^ (20 kpc) north-west of the CCG. 
The emission velocity exhibits a shear of ^ 250kms~^ along 
the same axis, possible due to rotation, with linewidths of ~ 
250 km s"^ FWHM. 

As for A 1664, the [NII]-i-Ha kinematic fits were used to search 
for associated emission in [SIIJAA6717, 6731. Unfortunately, the 
emission in [01] AA6300, 6363 is badly affected by night sky emis- 
sion line features and was not analysed. Line ratios as a function 
of local Ha flux are shown in Fig.[n] but the only firm conclusion 
that can be drawn is that the [NII]A6584/Ha is once again remark- 



ably constant with luminosity and consistent with the integrated slit 
value measured by Crawford et al. (1999). The [SII]A6731/6717 
ratio shows marginal evidence for an increase with Ha surface 
brightness, implying higher electron density as hinted at in A1664. 



A comparison with the CO emission is once again valu- 
able. The IRAM 30m single dish observation of Edge (2001) 
shows CO(l-O) from an inferred cool molecular mass of 1.8 di 
0.2 X 10^^ M©. The CO emission peaks at a velocity of 
— lOOkms"^ (relative to an assumed redshift of z — 0.2523), 
identical to that of the Ha emission peak. The CO linewidth of 
227 ± 38 km s"^ FWHM is again a close match to the Ha. Subse- 
quent Owen Valley Millimetre Array interferometry of the CO(l-O) 
emission constrains the angular extent of the CO to < 9" (36 kpc). 
It thus appears that, as in A1664, the CO and Ha emission share 
the same kinematics and morphology. 
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Figure 7. Fits to the Hck+[NII] complex in A 1664 for the regions labelled in Fig.|2| The two velocity components are shown in red and blue and the green line 
denotes their sum. Spikes adjacent to the B3 profile are sky emission feature residuals. 





Figure 9. Left: true-colour image of A1664 derived from archival VLT-FORSl V, R and I imaging. Centre: overlaid with contours of the red Ho; velocity 
component of Fig.|2| Right: overlaid with contours of the blue Ha velocity component. 



5 A2204 

Abell 2204 lies at ^ = 0.1514 and has an observed Ha slit lumi- 
nosity of 10^^ ergs~^ (uncorrected for reddening), although the 
inferred visible star formation rate is modest 1 M© yr~^ Craw- 
ford et al. 1999). The system is also a strong emitter of near-infrared 
H2 ro-vibrational emission (Edge et al. 2002), and the CO(l-O) 
emission from IRAM 30-m observations implies a molecular gas 
mass of 2.3 ± 0.6 x 10^° M© (Edge 2001). Chandra X-ray ob- 
servations reveal that the cluster has a disturbed core morphology, 
suggestive of a recent merger, with cold fronts at radii of ^ 28 kpc 



and 54.5 kpc, and the central radio source is also disturbed with 
three components within 10 arcsec, aligned roughly N-S (Sanders, 
Fabian & Taylor 2005). 

Results from the VIMOS data in Ha and the continuum are 
shown in Fig. [T2I along with an HST WFPC2 F606W (wide V- 
band) continuum image. The latter shows that the CCG is dusty and 
possibly interacting with galaxy #1, which is 15 kpc away in projec- 
tion at a relative velocity of 250 km (see Jenner 1974 and sec- 
tion 5.1). The Ha emission also has an irregular morphology with a 
number of filaments emanating from the nuclear regions, possibly 
coincident with the dust obscuration. The Ha velocity field has a 
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Figure 10. Continuum (66OO-68OOA, observed frame), Ho; intensity, velocity field and FWHM for A1835, as derived from single component fits to the 
Hq;-i-[NII] complex. The velocities are in units of kms~^ and relative to the nominal cluster redshift of z = 0.2523. As explained in the text, only part of 
the field of view of one of the two dithered pointings was used. 



tri-partite structure with three distinct clumps at velocities of ^- 
200, -50 and lOOkms"^ (relative to z = 0.1514). The velocity 
width of the Ha is in most locations ^ 200kms~^ FWHM, ex- 
cept in the nuclear region where it rises to 800 km . This in- 
creased velocity width is unlikely to be due to an AGN because the 
emission line ratios are uncharacteristic of AGN and show no spa- 
tial variation across the source. The Chandra X-ray data show no 
evidence for a nuclear point source. A similarly high velocity width 
in a spur to the north-east of the nucleus may be due to interac- 
tion between the emission line clouds and the radio source, whose 
core structure extends 5 arcsec from the nucleus in this direction. 
The high Ha line width may explain the apparent faintness of the 
CO(l-O) line (Edge 2001). VLT-ISAAC long-slit spectroscopy by 
Jaffe, Bremer & Baker (2005) of A2204 shows that the Paa and 
1-0 S(l) trace each other very closely, spatially and kinematically. 
Narrow-band Ha imaging by the same authors shows that the line 
emission drops off quite sharply, being characterised by a gaussian 



of FWHM ~ 30 kpc. Other clusters, in contrast, show a power-law 
emission extending out to radii of 40 kpc and beyond. 

As remarked above, the line ratios exhibit no spatial variations 
so are plotted as a function of local Ha flux in Fig. [T4| The con- 
stancy of the [NII]A6584/Ha and [OI]A6300/Ha is striking. As 
in the other targets, there is marginal evidence for an increase in 
[SII]A6717/Ha and [SII]A6731/6717 with increasing Ha surface 
brightness. 



5.1 Continuum sources in the field of A2204 

In Fig. ll5l we show spectra of the three continuum sources visible 
in the VIMOS field of view (Fig.[T3 alongside the CCG. The HST 
WFPC2 image of the same field shows that objects # 1 and 2 are 
galaxies and object #3 stellar. From the Na DA 5892.5 absorption 
blend we determined redshifts of z=0. 15096 (#1) and z=0. 13908 
(#2); the presence of z=0 Ha absorption confirms the stellarity of 
object #3. Galaxy #1 is thus close in projection (~ 15 kpc) and 
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Figure 12. Ho; intensity and 7000-7200 A (observed wavelength) continuum in A2204. The former is derived from a single component fit to the Hq;-i-[NII] 
complex. For comparison, we also shown an HST WFPC2 F606W image of the same region. In the latter, we identify three separate continuum objects distinct 
from the central cluster galaxy (CCG) and a z = 1.06 lensed background [Oil] emitter (object #4) (shown in a magnified portion of the HST frame in Fig.ll6i. 
all of which are discussed further in section 5. 



velocity (c^ 100 km s"^ ) to the CCG of A2204. It is thus possi- 
ble that the two are interacting. Galaxy #2 is twice as far away in 
projection and offset by -3200 km in velocity. 



a follow-up with a longer IFU observation at higher spatial reso- 
lution, possibly with adaptive optics due to the proximity of the 
moderately bright star #3 (R=15.3; B=16.5 mag). 



5.2 The gravitationally-lensed background galaxy 

The HST image of the cluster in Fig. fT2l exhibits a gravitationally 
lensed arc which terminates on a bright knot of emission which we 
refer to as object #4. The VIMOS data show [OII]AA3726, 3729 
emission in 1 1 fibres around this position which, with a redshift of 
z = 1.0604, we identify as the lensed background galaxy. Fig.fT6l 
shows the portion of the calibrated 2D spectral frame in which this 
emission was discovered, together with the reconstructed emission 
line image. The line is just at the edge of the 7600-7700A telluric 
absorption, for which we correct with the aid of the stellar con- 
tinuum spectrum of object #3. Line emission is present along the 
full length of the arc as defined by the HST image. The total [Oil] 
flux from it is estimated to be 7.2 x 10~^^ ergcm"^ , a fig- 
ure obtained by comparing our Ha measurments of the CCG with 
the slit flux in Crawford et al. (1999). The velocity of the emission 
varies by less than zb20kms~^ , and the line width is the range 
100-150 km s"^ FWHM. This lensed arc would be suitable for 



6 ZW 8193 

The core of Zw 8193 is complex in both space and velocity and 
initial optical spectroscopy of this cluster yielded discrepant red- 
shifts for the emission (z = 0.1829 ± 0.0001) and absorption line 
(z = 0.1725 ± 0.0001) components (Allen et al. 1992). Fig.[T7l 
shows that there are three galaxies within 2-3 arcsec of the CCG; 
UKIRT long-slit spectroscopy revealed Paa emission in and be- 
tween the CCG and its neighbour to the north, with an apparent 
velocity difference of 400 - 500 km s"^ (Edge et al. 2002). The 
benefit of higher spatial and spectral resolution OASIS IFU data 
now reveals the full complexity of the system. The Ha emission 
consists of three distinct 'blobs', one centred on the CCG, one 
on a galaxy to the north and one slightly offset from the latter 
to the west. The velocity ranges over zb250kms~^ (relative to 
z = 0.1825) and reconstructing the placement of the 1.2^^-wide 
N-S slit used for the UKIRT observations shows the consistency 
of the Paa and Ha velocity fields. The [NII]/Ha ratio varies from 
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Figure 13. Velocity field and FWHM (in kms~^ ) derived from fits to the Hq;+[NII] complex in A2204; the Ho; intensity map of Fig.[T2lis reproduced here 
for ease of comparison. The zero-point of the velocity field is taken to be the nominal cluster redshift of z = 0.1514. Line profiles for some individual fibres 
are also shown. 



0.9 on the CCG to 0.7 at the position of the northern blobs. The 
ratio [SII]A6717/Ha ~ 0.5 and is also relatively constant across 
the emission. The [SII]A6731 line is badly affected by a night- sky 
emission line. The extinction-corrected Ha slit luminosity quoted 
by Crawford et al. (1999) is 3.1 x 10^^ ergs"^ (converted to the 
cosmology used here). 

The CO(l-O) observations of Zw 8193 have yielded seemingly 
inconsistent results (Edge 2001). One observation revealed a broad 
plateau of emission from v = —100 to ^; = +500kms~^ but 
a subsequent observation failed to confirm it. As remarked by 
Edge (2001), the presence of a flat- spectrum radio source in 
Zw 8193 complicates the baseline subtraction. The quoted upper 
limit on the molecular gas mass is 4.3 x 10^° Mq. 

As in A1664, and to a lesser extent A2204, the Ha emission 
morphology and kinematics appear to be influenced by the close 
passage of a small cluster galaxy into the core. 



7 SUMMARY OF OBSERVATIONAL RESULTS 

Despite the diversity and complexity of these individual clusters 
and the small sample size, a number of generic features are appar- 
ent. They are: 

• A possible association between disturbed Ha emission and sec- 



ondary galaxies within kpc. In A1664 it appears plausible 

on both morphological and kinematic grounds that the Ha filament 
and complex nuclear structure in the CCG have been produced 
by the infall of a small cluster galaxy. This galaxy is now visible 
at the extremity of the filament. Since the mass of Ha-emitting 
ionized gas is relatively small (a few 10^ M©, as inferred from 
the Ha luminosity and the measured electron density), it is pos- 
sible that a significant fraction of it may have been stripped from 
the infalling galaxy. Our preferred interpretation, discussed in sec- 
tion 8, is that gas distribution has been disturbed by the infalling 
galaxy. In Zw8193 there is extended Ha emission associated with 
a galaxy which lies within 6 kpc (projected) of the CCG. In A2204, 
the irregular Ha kinematics may result from a recent interaction 
with a nearby galaxy (#1 in Fig. [T2l . or disturbance from the ir- 
regular radio source. In A1835, the Ha emission is elongated in 
the direction of a small companion galaxy 20 kpc away. Another 
well- studied CCG showing strong evidence for the influence of a 
nearby secondary galaxy on the Ha emission is RXJ 0820.9-1-0752 
at ^ = 0.110 (Bayer-Kim et al. 2002). 

• Uniformity of the ionization state. Despite the morphological 
and kinematic complexity of the Ha emission, the ionisation state 
of the gas is surprisingly uniform, as a function of position. Ha 
surface brightness and between clusters. The [NII]A6584/Ha and 
[OI]A6300/Ha ratios are particularly tight, but there may be hints 
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Figure 17. The Hq;+[NII] and [SII] line emission in Zw 8193 derived from the WHT/OASIS observations, compared with a K-band continuum image. North 
is up, east is to the left. 



of a weak increase of [SII]A6717/Ha and [SII]A6731/6717 with 
Ha surface brightness. These results suggest that the same mech- 
anism(s) excite all the optical line emission in these systems, and 
that they are quite independent of the processes which disturb the 
gas morphologically and kinematically on the kiloparsec scales we 
observe, i.e. an infalling galaxy or a weak radio source may stir up 
the gas, but they seem not to impact on its ionization state. For this 
reason, it seems unlikely that large-scale shocks are involved. 

The line ratios we measure are naturally consistent with 
the integrated slit values found by Crawford et al. (1999) and 
characteristic of high Ha luminosity (> lO^^ergs"^) cooling 
flow CCGs. The latter show strong AGN-like [NII]A6584/Ha 
and [01] A6300/Ha, but weak [SII] A6717/Ha and [OIII] A5007/H/3 
similar to starbursts. As mentioned in section 1, it was shown by 
Crawford et al. (1999) that the observed Ha luminosities of such 
systems can be powered by photoionisation by the 0-star popula- 
tions which explain their excess blue light. Typically 10^ — 10^ 
stars and larger populations of B, A and F stars are required in 
excess of the continuum expected from a cluster elliptical. How- 
ever, the ro- vibrational H2 :Paa line ratios of CCGs are much higher 
than typical starburst galaxies (Jaffe, Bremer & van der Werf 2001; 
Edge et al. 2002). This may require a distinct population of much 
denser clouds (Wilman et al. 2002) or excitation by a spectrum 
harder than that of 0-type stars (Jaffe, Bremer & Baker 2005). 
This may also explain why the optical emission line ratios differ 
from those of starburst galaxies. For the ensemble of line-emitting 
CCGs, Crawford et al. (1999) discovered a continuous trend of de- 
creasing [NII]A6584/Ha with increasing Ha luminosity, and simi- 
lar weaker trends in [OIII]A5007/H;^ and [SII]A6717/Ha. Our re- 
sults show that similar trends do not extend to local Ha surface 
brightness within the most luminous systems. 



Sparks et al. (1989) proposed that cool gas (e.g. from a galaxy 
merger) falling into the core of a galaxy cluster could be energised 
by thermal conduction from the X-ray-emitting ICM. This would 
lead to optical line emission from the cooler gas and a local en- 
hancement in X-ray emission. Such a model was used to account 
for the close correspondence between the X-ray and Ha filaments 
in the Virgo and Perseus clusters (Sparks et al. 2004). It is not 
known whether this mechanism can reproduce the observed emis- 
sion line ratios. 

• CO and Ha trace the same gas, at least kinematically and proba- 
bly also spatially. Two of the four clusters, A1664 and A1835, have 
robust detections of CO(l-O) emission in the survey of Edge (2001). 
In both these clusters the CO emission profile is in very good agree- 
ment with the Ha kinematics, in terms of its centroid velocity 
and line width, strongly suggesting that the CO and Ha emission 
arise from the same ensemble of clouds. Given this strong CO: Ha 
link, the weak (A2204) or ambiguous (Zw8193) CO detections in 
Edge (2001) most likely stem from a high Ha velocity width or 
zero-point offset (relative to the CCG redshift then assumed), re- 
spectively. In the aforementionned RXJ 0820.9-^0752 at z = 0.110 
(Bayer- Kim et al. 2002) the CO and extended Ha kinematics are 
also well-matched. This is in line with the correlation between Ha 
luminosity and molecular gas mass discovered by Edge (2001). In 
contrast, the ro-vibrational H2 emission and kinematics appear to 
be less strongly coupled to the Ha. In A1664, the H2 has a much 
lower velocity width than Paa and CO, although in A2204 the 
match is much closer. RXJ 0820.9-1-0752 has weak or non-existent 
H2 emission in the Edge et al. (2002) survey, despite its strong ex- 
tended CO and Ha emission. 

We note, however, that Jaffe, Baker & Bremer (2005) discov- 
ered that the Paa and hot H2 trace each other spatially and kine- 
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Figure 6. The 10 ks Chandra ACIS-S image of A1664 in the 0.5-5 keV en- 
ergy band with 1.5^^ gaussian smoothing (courtesy A.C. Fabian). The or- 
ange cross-hair denotes the position of the galaxy #1 (Fig.|2| close to the 
CCG. 
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Figure 5. lonisation diagnostic line ratios plotted as a function of Ho; flux 
for individual fibres in A 1664. Results are shown separately for the red and 
blue velocity components (where fitted), and also for fibres where a single 
velocity component was fitted. 



Figure 8. The Ho; and CO(l-O) emission line profiles integrated over com- 
parable areas of A1664. The CO(l-O) is derived from the IRAM 30m obser- 
vation by Edge (2001) with a 23.5^^ beam size (the solid and dashed lines 
show measurements with different correlators). The Ha profile was con- 
structed by summing over the whole source and using the multi-component 
gaussian fits to the Hq;-i-[NII] complex- the soHd shows all the emission, 
whilst the dotted and dashed line show the contributions from fibres re- 
quiring one and two velocity components, respectively. Velocities are with 
respect to the nominal cluster redshift of z = 0.1276. 
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Figure 11. Ionization diagnostic line ratios plotted as a function of Ho; flux 
for lenslets in A 1835. 
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Figure 14. Ionization diagnostic line ratios plotted as a function of Ho; flux 
for lenslets in A2204. 



matically out to radii beyond 20 kpc, but A2204 is the only cluster 
common to our two samples. 



8 INTERPRETATION: PHYSICAL PROCESSES IN THE 
DENSE MOLECULAR GAS RESERVOIR 

We now interpret the above findings, taking as our starting point the 
CO observations. As reviewed in section 1, these have revealed sub- 
stantial masses of cool molecular gas (lO^"^^'^ M©, at 20-40 K) 
in the cores of cooling flows, with a mass which correlates with 
the Ha luminosity (Edge 2001). Subsequent interferometry has 
shown that the CO emission is confined to scales < 20 kpc (Edge & 
Prayer 2003). The lack of double-peaked CO line profiles suggests 
that the molecular cloud distribution is quasi- spherical, rather than 
disk-like. Such a dense, compact molecular gas reservoir will most 
likely be disturbed by the infall of a small galaxy into the cluster 
core, as we believe is happening most clearly in A 1664, or by the 
onset of a nuclear radio source. We consider these possibilities in 
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Figure 15. Spectra of continuum objects #1-3 in the VIMOS field of view 
of A2204 (see Fig.IT2|. Redshifts were determined from the Na DA5892.5 
absorption blend (#1 and 2) and confirm their cluster membership. Object 
#3 is a star. These spectra have not been flux calibrated in either an absolute 
or relative sense. The features between 686O-688OA and 7600-7700A are 
telluric absorption. 



more detail below. First we consider some physical processes in 
the unperturbed molecular gas reservoir and their observable con- 
sequences. 

It is plausible that the CO and Ha emission arise within the 
same clouds, because they appear to share the same kinematics 
and are closely correlated in luminosity. The ionised surfaces of 
such clouds are likely to be in pressure equilibrium with the hot 
X-ray gas for which n ^ 0.1 cm~^, T ~ 10^ K; for optical 
emission line regions n ^ 100 cm~^ , T 10^ K, such that 
nT is the same for both phases. Pressure equality may also ex- 
tend to the cool molecular gas if n ^ 10^ cm~^ , although it is 
possible that the cloud cores are self-gravitating. The clouds may 
be warmed/ionised by a combination of young stars and X-ray 
emission/conduction from the intracluster medium, but the pre- 
cise details do not concern us here. The rovibrational H2 emis- 
sion may be a distinct, transiently-heated high-pressure component, 
since line ratio analysis points to thermal excitation in dense gas 
with n > 10^ cm"^ and T 2000 K (Jaffe, Bremer, van der 
Werf 2001; Wilman et al. 2002). Wilman et al. (2002) proposed 
that shocks between a population of low-density 200 cm~^ ) 
Ha-emitting clouds could be used to transiently create such high- 
density molecular gas, although an external stellar radiation field 
was still used to excite the observed emission. We now modify this 
model and investigate collisions amongst the denser CO clouds. 
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Figure 16. The z = 1.0604 gravitationally-lensed background galaxy in 
the field of A2204. (a) Part of the HST WFPC2 F606W frame showing 
the arc in the continuum; objects 2-4 are labelled as in Fig.[l2| (b) Recon- 
structed image in the [Oil] AA3726,3729 doublet using the VIMOS data; the 
surface brightness is in units of lO"-*^^ erg cm~^ s"-*^ arcsec"^. (c) A por- 
tion of the wavelength-calibrated 2D VIMOS spectral frame of A2204 in 
which the arc was discovered; boxes higlight the [OII]AA3726,3V29 emis- 
sion. 



8.1 Cloud-cloud collisions and external perturbations 

Consider that a mass M of CO-emitting molecular gas is confined 
within a radius Ro of the CCG, distributed within an ensemble of 
identical clouds of radius r, mass mc, space density J\f{lZ) and 
hydrogen number density n. Interferometry shows that the CO is 
centrally concentrated so we take N{Tl) oc Tl~^ , as in the inner 
part of an NFW profile (Navarro, Frenk & White 1997). For a typi- 
cal relative cloud velocity of v, the rate of cloud collisions per unit 
volume is J\f{lZY crv, where a is the collision cross-section, which 
we approximate with the geometrical cross-section a — nr^. It 
then follows that the rate at which mass is processed through such 
collisions is: 



McOL 



SM^v 



(1) 



For V 



150kms"\ Ro = lOkpc, 



: Ipc and n = 

' , eqn. (1) evaluates to Mcoll ^1-100 M© yr~^ for 
M — 10^°~^^ M© (mp is the proton mass). Such supersonic col- 
lisions will shock the clouds to temperatures of several 10^ K. 
For the assumed densities the cooling time from these tempera- 
tures is very short, 1 month, so the gas will quickly return 
to the cool phase, i.e. collisions will merely recycle the cool gas 
and not deplete it. It is, however, possible the shock compression 



will act as a trigger for star formation. They may also create the 
dense molecular gas and a locally strong, hard (e.g. black body 
at T > 10^ K) radiation field required for the production of the 
rovibrational H2 emission. The EUV luminosity produced by such 
shocks is LcoLL = ^Mcollv'^, or: 



An'^R^nrUpr ' 



(2) 



For the same parameters as before, this amounts to Lcoll = 
7 X io^^~^^ ergs~^ . Although this is several orders of magni- 
tude below the inferred luminosity of the young stellar component 
in these systems, the cubic dependence on v and Ro and uncer- 
tainty over the actual values of r and n means that such shocks 
could make a non-negligible contribution to the excitation of the 
line emission in CCGs, at least locally if the gas clouds are not 
uniformly distributed. We note, however, that limits on the strength 
of the [OIII]A4363 emission line imply that shocks are unlikely to 
play a strong role in the production of the optical line emission in 
Ha-luminous CCGs (e.g. Voit & Donahue 1997). The mechanism 
may, however, be a relatively more important source of excitation 
in lower-Ha luminosity systems. Such shocks will steadily drain 
energy from the molecular gas reservoir, forcing the distribution to 
smaller radii and thus increasing the collision rate still further. 

We now consider the consequences of disturbing this reser- 
voir with the infall of a small galaxy, as appears to be happening 
most dramatically in A 1664. From the constant excitation state of 
the ionised gas, it seems that this does not cause any large scale 
shocks or cloud destruction. Most likely, the cloud distribution will 
be dynamically perturbed by the infalling galaxy, producing the ob- 
served streams and filaments. Such interactions are likely to en- 
hance the star-formation rate, either by increasing the cloud-cloud 
collision rate (as above), or via some other means. With reference 
to the relationship between [NII]/Ha and Ha luminosity discov- 
ered by Crawford et al. (1999), the observed uniformity of the line 
ratios we find in our sample may be due to the fact that, at high 
Ha luminosity, all parts of the CCG are saturated at the lowest 
[NII]/Ha values characteristic of massive star formation. In less 
Ha-luminous CCGs, we would thus predict more variations in the 
line ratio, with [NII]/Ha increasing away from the less-abundant 
star-forming clumps. 

The possibility that galaxy interactions may trigger star- 
formation in these molecular gas-rich environments highlights the 
similarities between these CCGs and luminous infrared galaxies, 
as discussed by Edge (2001). The Spitzer observations of Egami et 
al. (2006) further strengthen the connection, revealing large thermal 
dust luminosities Lir >lO^^L0in two of the most Ha-luminous 
CCGs (Zw3146 and A1835). 



9 CONCLUSIONS 

The results presented in this paper have shed new light on the kine- 
matic and morphological complexity of the emission line gas in 
the cores of these Ha-luminous clusters. Our principal findings are 
that: • In many cases the Ha emission appears to have been dis- 
turbed by interaction with a small companion galaxy to the central 
galaxy. • The CO and Ha emission share the same kinematics and 
by implication trace each other closely. • The emission line ratios 
of the gas are uniform, showing no strong variations in [NII]/Hq;, 
[SII]/Ha or [OI]/Ha, as a function of position or Ha surface bright- 
ness, suggesting that its excitation state is independent of the pro- 
cesses which disturb it on the observed kiloparsec scales. 
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Our interpretation is that the Ha and CO emission are pro- 
duced in a population of molecular clouds warmed by a starburst, 
which has itself been triggered by the passage of a secondary 
galaxy through the gas-rich cluster core. In these most luminous 
Ha CCGs, this galaxy-wide starburst leads to the saturation of 
the optical emission line ratios at a uniform level characteristic 
of massive star formation. In the absence of such triggering, the 
cooled gas reservoir will be in a more quiescent state with a much 
lower level of Ha and CO emission; optical emission line ratios 
will be governed by other processes, e.g. excitation through cloud- 
cloud collisions (as described above) or mixing layers (Crawford 
& Fabian 1992), to mention just two possibilities. Starting from 
this low Ha luminosity quiescent state, adding some star-forming 
clumps will increase the overall Ha and CO luminosity and give 
rise to spatial variations in the optical line ratios. At the highest Ha 
luminosities, star-formation will be widespread, with little spatial 
variation in optical emission line ratios. 

To test this interpretation, future IFU observations should be 
carried out on CCGs covering the [Oil] A3727 and [OIII]-hH;^ emis- 
sion lines and the sub-4000A continuum in order to probe current 
star formation, and with near-infrared IFUs to study the H2 and 
Paa emission. Observations should also be performed on lower Ha 
luminosity CCGs in order to examine whether spatial variations in 
optical emission line ratios are indeed larger, due to less widespread 
star formation. On a longer timescale, mm-interferometers such as 
ALMA will reveal the spatial distribution of the CO emission (and 
dust) for a full comparison with the Ha. Our findings also high- 
light the need for numerical modelling of these dense CO reser- 
voirs, with and without external perturbations. 
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